The ultimate state of thermal convection 
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The ultimate regime of thermal convection, the so called 
Kraichnan regime (R. H. Kraichnan, Phys. Fluids 5, 1374 
(1962)), hitherto has been elusive. Here, numerical evidence 
for that regime is presented by performing simulations of the 
bulk of turbulence only, eliminating the thermal and kinetic 
boundary layers and replacing them by periodic boundary 
conditions. 



Thermally driven turbulence is of tremendous impor- 
tance in oceanography, geophysics, meteorology, astro- 
physics, or process technology. If the flow is confined to a 
box heated from below and cooled from above, the setup 
is called Rayleigh-Benard convection. It is one of the 
classical problems in fluid dynamics [|],^. The control 
parameter are the dimensionless temperature difference 
between bottom and top (called Rayleigh number Ra), 
the ratio between kinematic viscosity v and heat diffusiv- 
ity K (called Prandtl number Pr), and the aspect ratio 
of the cell. The system answers with some heat transfer 
from the bottom to the top (called Nusselt number Nu 
in dimensionless form) and the degree of turbulence in 
the cell (expressed in terms of the Reynolds number Re). 

All classical experiments on this system show (effec- 
tive) power law dependences with Nu ~ RaP "^^ up to 
Nu ~ RaP'^^. However, for very large Ra, Kraichnan 
1^ has predicted an ultimate scaling law Nu ^ Ra}/"^. 
In this ultimate state the heat flux and the turbulent 
velocity are expected to be independent of the viscos- 
ity and heat diffusivity. Indeed, for a system with in- 
finite aspect ratio a mathematically strict upper bound 
Nu < 0.167i?ai/2 _ I could be found [|. 

Hitherto, this ultimate state of thermal convection re- 
mained elusive inspite of tremendous efforts to find 
it. Niemela et al. |^ measured up to Ra « 10^* in he- 
lium gas close to the critical point (Pr « 1), still finding 
an effective power law Nu ~ Ra^ "^^. Glazier et al. 
measured up to Ra « 10"'^^ in mercury {Pr — 0.025), 
finding Nu ^ Ra'^^^^. The only evidence for a transi- 
tion (at Ra sa 10^^) towards the ultimate regime has 
been claimed by Chavanne et al. [^|^, but Niemela and 
Sreenivasan argue that those data would be consis- 
tent with Nu ~ Ra^/^. In summary, it is therefore not 
clear 

• whether in today's experiments the Rayleigh num- 
ber is only still too small for the ultimate state, 
or 



• whether it does not exist at all. 

However, it is crucial to know which of these two al- 
ternative is the correct one: 

• for practical reasons in view of above mentioned ap- 
plications of thermal convection in geophysics, in- 
dustry etc, as upscaling from lab-scale models does 
not work if there is a transition to a new state of 
turbulence; 

• for fundamental reasons, as is a "central dogma 
in turbulence" ^ that its effects become indepen- 
dent of viscosity and diffusivity for large enough 
Reynolds (or here Rayleigh) numbers pi[ . 

Indeed, also within Grossmann and Lohse's unifying the- 
ory on thermal convection ||l^, |l3| the Kraichnan or ulti- 
mate regime follows, once the thermal and kinetic bound- 
ary layers either break down due to their expected in- 
stability at large Ra, or hardly contribute to the global 
kinetic and thermal dissipation |jl^. According to ref. 
p3[ the breakdown of the laminar kinetic boundary layer 
should happen at Ra « lO^'' for Pr — 1 and at Ra « 10^^ 
for Pr — 0.025 for an aspect ratio F = 1 cell; Niemela 
and Sreenivasan [|l^ find similar values. Then the total 
kinetic dissipation rate e„ equals the bulk kinetic dis- 
sipation eu,buik and the total thermal dissipation rate 
eg equals the bulk thermal dissipation ee,buik and the 
Kraichnan or ultimate regime with 



and 



Nu - Ra'^I^Pr^l^ 



Re - Ra'^/^Pr"'^/^ 



(1) 



(2) 



results (ig. 

From the above it follows that an artificial destruc- 
tion of the boundary layers should enhance the ultimate 
regime. Indeed, by performing Rayleigh-Benard convec- 
tion experiments in a cell with rough top and bottom 
walls Roche et al. find evidence for the onset of the 
Kraichnan regime. 

Here, we follow another approach: We perform numer- 
ical simulation o/ the bulk of turbulence only, eliminating 
the thermal and kinetic boundary layers and replacing 
them by periodic boundary conditions [ p6| . Then the 
Kraichnan regime should follow immediately. The ther- 
mal bulk turbulence is forced by a mean temperature 
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FIG. 1. (a) Nu (squares, lower curve) and Re (circles, up- 
per curve) as a function of Ra from our simulation. The 
error bars are based on the variances of measured quanti- 
ties which have been obtained by splitting the total time se- 
ries in smaller pieces. Note that in this way we could also 
check that our simulations are statistically stationary. The 
straight lines are linear regressions giving Nu ~ ija"'^^*""® 
and i?e ~ i?a0-55±0 02 

(b) Same as in (a), but now compensated by Ra^^^. 



profile —zA/L, where L is the periodicity length and A 
the temperature difference on that length scale. The tem- 
perature fluctuations 6 then obey the advection equation 



dt9 + UidiO = ndidiO + u^- 



(3) 



The velocity field Ui{x, t) obeys the standard Boussinesq 
equation, 



(4) 



(3 is the thermal expansion coefficient, g gravity, and p 
the pressure. Ra and Pr are defined as usual as Ra — 
l3gAL^/{i/K) and Pr = respectively. The Prandtl 
number is 1 in our simulations presented here. Also Nu 
is defined as usual as Nu — {u^9) L/{kA) — (0)^ L/A. 
The Reynolds number is defined with the rms mean ve- 
locity fluctuation u' , i.e.. Re = u'L/v. A similar simu- 
lation has earlier been performed by Borue and Orszag 
[p^ . The focus of that work was on scaling properties 
and spectra and thus hyperviscosity has been employed. 



From eqs. (y) and <^ one can derive by volume averag- 
ing and assuming statistically stationarity the exact p7| ] 
relations eg — kL^^A^Nu and e„ = u'^ N uRaPr . 
As in refs. we estimate Cufiuik ~ u'^/i and ee^uik ~ 
w'A^/L and the Kraichnan regime (0) - (|) results. Note 
that L is the relevant length scale for these estimates as 
for the Bolgiano length scale it holds > L. 

This Kraichnan scaling is indeed seen in our numerical 
simulations of eqs. (|) and (|). We run four different Ra 
numbers Ra = 8.64 • 10^ (for 160 large eddy turnovers), 
Ra = 2.21 ■ 10^ (for 60), Ra = 4.51 ■ 10^ (for 180), and 
Ra = 1.38-10^ (for 334), and average over space and time 
to obtain Nu and Re The i?a-scaling is consistent 
with eqs. ([|) - (§), see figure]^. To our knowledge it is the 
first realization of the ultimate regime of thermal convec- 
tion. It is remarkable that it is realized in spite of the 
relatively low Ra numbers of the simulations. The reason 
is that the boundary layers have been eliminated and the 
simulations focus on the bulk. Due to the lack of bound- 
ary layers and the more efficient driving, the Reynolds 
numbers achieved here are much larger than they would 
be in the standard Rayleigh-Benard case. Vice versa, the 
required Rayleigh numbers to achieve the same degree 
of turbulence for the standard Rayleigh-Benard case as 
here are much higher. E.g., to achieve Re — 4996 (here 
for Ra = 1.38 ■ 10'') one needs a Rayleigh number of 
Ra ^ 7.2 ■ 10^ in the standard case f^. 

In conclusion, from above alternative the first one is 
favorable: Once Ra is so large that the laminar kinetic 
boundary layers breaks down, the ultimate regime of 
thermal convection, which has been so elusive in exper- 
iment, should finally show up. Obviously, the ultimate 
proof for the existence of this regime can only come from 
experiment. 
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